epl draft 



Origin of intermittency in wave turbulence 



E. Falcon 1 W, S. G. Roux 2 and C. Laroche 1 



1 Laboratoire Matiere et Systemes Complexes (MSC), Universite Paris Diderot, CNRS - UMR 7057 
10 rue A. Domon & L. Duquet, 75 013 Paris, France, EU 
. 2 Laboratoire de Physique, Ecole Normale Superieure de Lyon, CNRS - UMR 5672, 69 007 Lyon, France, 



EU 



o. 

(NT 



a: 

• »— i ■ 

CZ2 
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Abstract. - Using standard signal processing tools, we experimentally report that intermittency 
of wave turbulence on the surface of a fluid occurs even when two typical large-scale coherent 
structures (gravity wave breakings and bursts of capillary waves on steep gravity waves) are not 
taken into account. We also show that intermittency depends on the power injected into the 
waves. The dependence of the power-law exponent of the gravity-wave spectrum on the forcing 
amplitude cannot be also ascribed to these coherent structures. Statistics of these both events are 
studied. 



^V Understanding the origin of intermittency is a challeng- 
f^ ing problem in varied domains involving turbulent flows. 
i '•~ faitermittency is the occurrence of bursts of intense mo- 
tion within more quiescent fluid flow [TJ[5]. This leads 
,- ^ to strong deviations from Gaussian statistics that become 

^ larger and larger when considering fluctuations at smaller 
^ I and smaller scales. In three-dimensional hydrodynamic 
„—*- turbulence, the origin of these deviations has been as- 
C^) cribed to the formation of coherent structures (strong vor- 
iy-!tices) since the 50's pQ. However, the physical mechanism 
("***) of intermittency is still an open question [3] . Intermit- 
^^tency has also been observed in granular systems [4], in 

"^ magnetohydrodynamic turbulence in geophysics J5] or in 

£> the solar wind [5] , and in systems involving transport by a 
? *~j turbulent flow [7] . A recent observation of intermittency 
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has been reported in wave turbulence [5], a system that 



C^ strongly differs from high Reynolds number hydrodynamic 
turbulence. It could thus motivate explanations of inter- 
mittency different than the ones considering the dynamics 
of the Navier-Stokes equation. 

The aim of this Letter is to investigate if some coherent 
structures are responsible of intermittency in wave turbu- 
lence. In the case of wave turbulence on a surface of a 
fluid, coherent structures such as bursts of capillary waves 
on steep gravity waves |Hj and wave breakings |10| are well- 
known phenomena, these latters being recently taken into 
account in numerical simulations |11) . Wave breakings 
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also occur in plasma waves, internal waves, and Rossby 
waves in geophysics. It has been suggested that inter- 
mittency in wave turbulence may be connected to wave 
structures (such as cusps, whitecaps or wave breakings) 
thus motivating theoretical [T2] and numerical [T3J works. 
Here, we show experimentally that intermittency does not 
come from wave breakings and capillary bursts on gravity 
waves. Using standard signal processing tools, one finds 
criteria to detect such structures that allow us to study 
their statistics and their possible role in the origin of in- 
termittency. We also show that intermittency depends on 
the power injected into the waves. The frequency- power 
law of the gravity-wave spectrum is known to depend on 
the forcing parameters [T4] • We show that this dependence 
is not related to these coherent structures. 

The experimental setup has been described previously 
[14J. It consists of a square vessel, 20 x 20 cm 2 filled 
with mercury up to a height of 2.6 cm. Similar results 
are found with water. Surface waves are generated by 
the horizontal motion of a rectangular plunging plastic 
wave maker driven by an electromagnetic exciter. This 
vibration exciter is driven with a random forcing within 
a narrow low- frequency range (typically 0.1 to 5 Hz), and 
a rms voltage amplitude ajj from 0.1 to 0.8 V leading to 
wave mean steepnesses (ratio of crest-to-trough amplitude 
to its duration) from 1 up to 4 cm/s. The rms value ay of 
the velocity fluctuations of the wave maker is proportional 
to o-jj. The mean power injected (J) into the fluid scales 
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Fig. 1: (color online) Temporal evolution of the surface wave 
amplitude (top) and the local slope of the wave (bottom) for a 
strong forcing (ou = 0.8 V) during 5 s. 



as (I) ~ <Ty ~ Oy [T3]. £ = (I)/(pA) is the mean energy 
flux where p is the fluid density and A the immersed area 
of the wave maker. The surface wave amplitude, rj(t), is 
measured at a given location of the surface by a capacitive 
wire gauge [T3]. r)(t) is low-pass filtered at 1 kHz and 
recorded with a 4 kHz sampling rate during 3000 s, leading 
to N = 1.2 x 10 7 points. This signal is cut in 10 files of 
300 s. Statistical properties of each file are then computed, 
the rms value of the computed quantity giving its error 
bar. This also allows us to check the signal stationarity. 

For a weak forcing, the wave amplitude r](t) is found to 
fluctuate around a zero mean value in a roughly Gaussian 
way, as well as the local slope of the surface waves 5 n (t) 
computed from the differential of r)(t). For a higher forcing 
amplitude, a typical temporal recording of the wave ampli- 
tude is displayed in Fig.[TJ rj(t) much more fluctuates with 
more probable high crest waves than deep trough waves. 
This comes from nonlinear effects due to the strong steep- 
ness of the waves. The corresponding local slope 6 v (t) is 
also strongly erratic (see Fig. [T]), and two typical events 
can be observed: short peaks of very high amplitudes, and 
trains of oscillations of much smaller amplitudes (see Fig. 
[T]and below). Both events occur randomly, and are always 
found close to the maximum of the local slope of the wave. 

Figure [2] shows such typical events that are detected on 
the wave amplitude signal once the forcing is high enough. 
Figure^ shows the first kind of event: a burst of high- 
frequency capillary waves on a gravity wave. This is a 
well-known phenomenon occurring when the gravity wave 
is steep enough [5]. Indeed, when the gravity wave ampli- 
tude increases, the local curvature at the crest increases 
rapidly leading to strong surface-tension effects close to 
the crest. A train of capillary waves then is emitted prop- 
agating down the gravity wave front face as predicted the- 
oretically [16] , One can also observe in Fig. [5^ that the 
wave is much steeper in the front than in the rear. The 
capillary wave frequency (obtained from the period be- 
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Fig. 2: (color online) Typical events occurring near the crest of 
a steep gravity wave: a) Capillary waves generation (ajj = 0.4 
V). b) Wave breaking (au = 0.8 V). Both events occur on the 
forward face of the gravity wave (wave front is left-hand side) . 
Axis of the wave amplitude (resp. local slope) is on the left- 
hand side (resp. right-hand side). Upper (resp. lower) curve 
corresponds to the wave amplitude (resp. wave slope). 



tween two successive peaks in the slope trace) is found to 
increase with the index of number of capillary waves from 
the gravity wave crest as already shown experimentally [3] 
and predicted theoretically [175]. The typical frequency of 
the carrier gravity wave is of the order of 5 Hz, whereas 
the capillary wave ones are in the range 80 — 250 Hz. 

Figure [5p shows the second type of event observed at 
high enough forcing. Sharp peaks occur on the wave-slope 
signal, corresponding to the early stage of a wave break- 
ing: a "bulge" is formed on the forward face near the crest 
preceded by small-amplitude capillary waves. It is well- 
known that as the wave steepens, the amplitude of the 
bulge increases [TU]. The leading edge of the bulge (also 
called the "toe") marks the formation of a train of small- 
amplitude capillary waves. These capillary waves gener- 
ally grow rapidly with time leading to the breaking of the 
wave |10j. 

As shown in Fig. [51 a wave breaking has a larger slope 
amplitude than the one of a burst of ripples generated on 
the gravity wave. Consequently, one can find a criterion to 
detect wave-breaking events in order to study their statis- 
tics. To wit, the local acceleration of the surface waves is 
computed from the second differential of rj(t). The proba- 
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Fig. 3: (Color online) PDF of the normalized wave acceler- 
ation, acc/oacc, for different forcings ou = 0.2 (V), 0.3 (x), 
0.4 (o), 0.5 (o) and 0.8 (□) V (from lower to upper curves). 
Gaussian with zero mean and unit standard deviation (black 
dashed line). Wave breaking onset is ±4cr acc [red (light gray) 
dot-dashed lines]. Inset: a acc as a function of the forcing. 



bility density function (PDF) of the acceleration, ace, nor- 
malized to its rms value, cr acc , is shown in Fig.[3]for differ- 
ent forcing amplitudes. The core of the PDF is Gaussian 
and independent of the forcing up to a critical wave accel- 
eration of ±4cr acc . Above this value, the PDF tails depend 
on the forcing: the larger the forcing, the more probable 
the rare events are. This critical acceleration is the onset 
of wave breakings. Indeed, as directly observed on the ac- 
celeration signal, wave breakings occur when ace > 4er acc . 
One can thus remove from the acceleration statistics the 
wave-breaking events, i.e. a set of short signal durations 
St = tf — ti where the absolute value of the acceleration 
becomes larger at time U (resp. lower at tf) than this 
threshold (typically St ~ 100 ms). When removing these 
wave-breaking events, the PDF of the filtered acceleration 
is then found to be almost Gaussian whatever the forcing. 
Note that this should not be confused with a lack of inter- 
mittency (see below). The inset of Fig. [3] shows that a acc 
increases with the forcing as expected. 

Statistics of wave-breakings is then performed using the 
above detection criterion. Figure 5J: shows the number of 
wave breaking detected as a function of the forcing. For 
the smallest forcing, no wave breaking occurs on the fluid 
surface. When the forcing is increased, the number of 
wave breakings increases. Note that the PDF of time lag 
between two consecutive wave-breaking events is found to 
exponentially decrease as expected for Poissonian statis- 
tics (not shown here). The PDF of a number of events 
occurring in a fixed period of time is also found to follow 
a Poisson law with an occurrence rate of events increasing 
with the forcing. 

Let us now focus on the detection of capillary bursts on 
gravity waves. A time-frequency analysis based on wavelet 



transforms is a useful method to analyse signals with mul- 
tiple time- varying frequencies [17J. It provides temporal 
and spectral information simultaneously and is thus well 
adapted to detect capillary burst events by thresholding 
the energy (i.e., the wavelet coefficients squared) contains 
in the bandwidth 50 - 250 Hz. One thus obtains a set 
of signal durations St = tf — f, where the wavelets coeffi- 
cients are above the threshold. St is of the order of 40 to 
80 ms from the smallest to the strongest forcing. Figure 
[4J: shows the number of capillary bursts detected by this 
process as a function of the forcing. When the forcing is 
increased, the number of capillary bursts is found to in- 
crease. We also find that the statistics of capillary bursts 
does not follow a Poissonian distribution. 

Wave breaking and capillary burst events are two dif- 
ferent coherent structures that can independently be de- 
tected within the signal 77(f) (see above). One can 
thus probe their respective role in the intermittency phe- 
nomenon in wave turbulence. The intermittent proper- 
ties of a stochastic stationary signal are generally tested 
by computing the structure functions using the first-order 
differences of the signal 8t](t) = r)(t + r) — 77(f). However, 
when the signal has a steep power spectrum, E v (f) ~ f~ n 
with n > 3, high-order difference statistics is required 
[5]. For gravity wave turbulence, the theoretical expo- 
nent of the power spectrum of wave amplitude is n = 4 
|18| . whereas experimentally it is found to depend on the 
forcing with n > 4 |14| , the origin of the discrepancy being 
an open problem. We found that statistical convergence of 
the structure functions is reached when using the fourth- 
order (or higher) difference statistics. This is due to our 
locally multi-derivable signal [T§] . The fourth-order differ- 
ences of the signal At^t) = 77(f + 2r) — 4?y(f + t) + 677(f) — 
Ar)(t — t) + r/(t — 2t), are thus computed in the following. 

The PDFs of Arj t (t) normalized to their rms values <7atj 
are computed, for different time lags 20 < t < 130 ms, ei- 
ther from the whole signal 77(f) (see Fig. |4^i) or from the 
signal where both wave breaking and capillary burst events 
are removed (see Fig.[2J)). In both cases, a shape deforma- 
tion of the PDFs is observed with r. The PDF is roughly 
Gaussian at large r, and its shape changes continuously 
when t is decreased. This is a direct signature of inter- 
mittency {20]. Since this latter is observed in both cases 
(Fig. [3Ji anddJ}), it clearly means that wave breaking and 
capillary burst events are not the origin of intermittency. 

The structure functions are defined by 



1 N 
S p (r) ^ (\A Vt (r)n ^ -J2\^lt(rW 



(1) 



where N is the total number of points in the signal. S p (t) 
are shown in Fig. [5] for a fixed forcing. For 15 < r < 100 
ms, all the structure functions of order p (from 1 to 6) are 
found to be power laws of t: S p (t) ~ r^ p where £ p is an 
increasing function of the order p. The exponents £ p of 
the structure functions are then plotted in the main Fig. 
[5] as a function of p for different forcings. £ p is found to be 
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Fig. 4: (Color online) PDF of normalized increments A-q/aAri 
for different time lags r = 20, 23, 27, 32, 38, 51, 70, 132 ms 
(from top to bottom) computed from: a) whole signal, b) sig- 
nal without wave breakings and capillary burst events. Gaus- 
sian fit (dashed line). Curves have been shifted for clarity. 
ou = 0.3 V. c) Number of detected events vs. the forcing. 



a nonlinear function of p such that £ p = C\p — ^rp 2 where 
the values of c\ and C2 are found to both depend on the 
forcing (see o-symbols in top and bottom insets of Fig. [5]). 
The nonlinearity of £ p (c2 ^ 0) is a second signature of 
intermittency [20J. The so-called intermittency coefficient 
C2 is found to increase from 0.2 up to 0.4 when the forc- 
ing is increased whereas ci is found to decrease from 2.8 to 
2.2. Intermittency is observed here over almost one decade 
in time (15 < r < 100 ms), corresponding to frequencies 
5 < l/(2r) < 33 Hz related to gravity wave turbulence 
regime. Indeed, as observed on the power spectrum of the 
wave amplitudes (not shown here), this upper boundary 
value is the crossover frequency between gravity and cap- 
illary wave turbulence regimes for the driving frequency 
bandwidth used. This crossover is known to depend on 
the forcing parameters [T3], and is thus slightly increased 
with respect to the one between linear gravity and capil- 

/ 3 \l/4 

lary waves — =- I S-E J (~ 17 Hz for mercury). 

From the definition of the structure functions of Eq. ([1]), 
one can remove the coefficients fS.r) t (r) obtained at times t 
in a given neighborough of a wave-breaking or a capillary- 
burst event, i.e. for t e [t t -T, t f +T}. We choose T = 160 
ms to ensure that all kept coefficients are not polluted by 
the event. Moreover, T has to be larger than the max- 
imum time scale (100 ms) used for the estimation of £ p . 
The drawback of this method is that a lot of statistics is 
removed just for a single event. Using this signal process- 
ing protocol, one can thus performed a structure function 
analysis by removing either all the wave-breaking events 
detected (D-symbols in the insets of Fig. [5]), or by remov- 
ing both the wave breaking and capillary burst events (*- 




Fig. 5: (Color online) Structure functions S p (t) of the fourth- 
order differences of the wave amplitude, At/, as functions of 
the time lag r, for 1 < p < 6 (as labeled). Solid line: Power 
law fits, S p ~ r 4p , where the slopes £ p depend on the order p 
(see Fig. [§}. Curves have been shifted for clarity, ou = 0.3 V. 



symbols in the insets) from the statistics in Eq. (JIJ . When 
wave breakings are removed, the intermittency coefficient 
C2 is found to slightly decreases for a fixed forcing, but 
C2 is still non zero whatever the forcing (see D-symbols 
in the bottom inset). This means that intermittency is 
still observed, and thus cannot be ascribed to the wave 
breakings. For a fixed forcing, when both wave breaking 
and capillary burst events are removed, the Ci strongly 
decreases (^-symbols), but it is still above 0.1, that is, one 
order of magnitude larger than the typical values found 
in usual hydrodynamic turbulence [20j . Capillary bursts 
thus enhance intermittency but are not its primary ori- 
gin. We have checked that these results do not depend 
on the signal processing used. When changing T strongly 
(from 160 to 650 ms), intermittency is still observed (c2 
decreases slightly from 0.13 to 0.11 for the lowest forcing). 
This means that although from 15% to 40% of the original 
signal is removed, intermittency is very robust and could 
not be ascribed to possible missing events of the signal 
processing. To sum up, both PDF and structure function 
analyses lead to coherent results showing that intermit- 
tency occurs even when two typical coherent structures 
(wave breakings and capillary bursts) are not taken into 
account. 

Finally, let us now focus on the frequency exponent of 
the gravity wave spectrum that is related to the second 
order structure function by n = £2 + 1 = 2(ci — C2) + 1. 
From the values of ci and C2, our measurements show a 
decrease of n from 6.2 to 4.6 for an increasing forcing. 
When removing the coherent structures from the signal, 
the value of c\ is found to be independent of the pres- 
ence of these both events (see top inset of Fig. [5]), and 
the exponent n of the gravity spectrum, estimated from 
n = 2(ci — C2) + 1, is still found to be a decreasing func- 
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Fig. 6: (Color online) Exponents £ p of the structure functions 
as a function of p for different forcings: cru — (V) 0.2 , (x) 0.3, 
(o) 0.4, (o) 0.5 and (□) 0.8 V (from upper to lower curves). 
Solid lines are best fits £ p = c\p — ^-p 2 . £ p are computed 
from the fourth-order differences with 15 < r < fOO ms (e.g., 
see Fig. [5j). Top and bottom insets: evolution of c\ and ci 
with the forcing: (o) whole signal, (□) signal without wave 
breakings, and (*) signal without wave breakings and capillary 
bursts (for the three lowest forcings). 



tion of the forcing amplitude. Coherent structures are 
thus not the origin of the forcing-dependent frequency- 
power law spectrum. Such a forcing-dependent exponent 
is coherent with previous direct observations of the gravity 
wave spectrum [14,21 . This departure from the theoreti- 
cal spectrum of gravity wave E* heo (f) ~ £ 1 ^ 3 gf~ 4 (with g 
the acceleration of the gravity and e the mean energy flux) 
[18] could be related to finite size effect of the container 
and/or to strong nonlinear effect of high wave steepness in 
experiments |141I21| . Indeed, weak turbulence theory as- 
sumes infinite basin and weak nonlinearity |18| . For this 
weak regime, a simple dimensional analysis thus leads to 
5* /leo (r) ~ e p/ 6 gP/ 2 r 3 P/ 2 . Here, one finds experimentally, 
for small p, S p {t) - e ap r Cl P where a = 0.4 ± 0.05 and a 
depends on the forcing (see inset of Fig. [5]). 

In conclusion, we have shown that intermittency in wave 
turbulence persists when two typical large-scale coherent 
structures (wave breakings and capillary bursts on steep 
gravity waves) are removed from the wave amplitude sig- 
nal. The power-law exponent of the gravity-wave spec- 
trum is known to be non-universal in laboratory experi- 
ments |141I21| . Here, we show that this exponent depen- 
dence on the forcing parameters cannot be also ascribed to 
these coherent structures. The origin of the intermittency 
phenomenon in wave turbulence is still an open issue. It 
could be ascribed to the large fluctuations of the energy 
flux [T5] or to other possible wave structures. 
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